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PREFACE 


This  report  describes  the  auroral  borealis  research  effort  being  per- 
formed by  M&S  Computing,  Inc.  for  the  Army  Ballistic  Missile  Defense 
Command  in  Huntsville,  Alabama.  The  primary  objective  of  the  study  is  to 
advance  the  understanding  of  the  auroral  phenomenon,  especially  of  its  inter- 
actions with  radar.  The  primary  data  gathering  instrument  is  the  Safeguard 
Perimeter  Acquisition  Radar  in  North  Dakota.  This  radar  has  been  used  to 
collect  large  quantities  of  high  resolution  auroral  backscatter  data  with  the 
simultaneous  tracking  of  a number  of  selected  satellites. 

The  data  presented  here  represents  the  initial  results  of  the  reduction 
of  data  collected  on  the  two  dates  of  September  17-18  and  September  26-27, 
1975.  The  report  provides  a detailed  description  of  the  data  collection  process 
and  the  data  reduction  techniques,  as  well  as  future  auroral  study  efforts. 
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INTRODUCTION 


The  purpose  of  this  report  is  to  describe  the  auroral  experiments 
performed  by  BMDSCOM  at  the  SAFEGUARD  PAR  Site  in  North  Dakota 
during  the  fall  of  1975.  The  objective  was  to  gather  a large  quantity  of  high- 
resolution  auroral  data  for  in-depth  analysis.  The  region  of  observation  was 
very  large  and  the  data  was  gathered  continuously  from  near  the  onset  of 
auroral  activity  until  it  ended.  In  addition  to  the  radar  data  tapes  obtained  at 
the  PAR  Site,  a collection  of  video  tapes,  charts,  graphs,  and  associated 
data  were  obtained  from  several  outside  agencies.  This  data  includes  geo- 
magnetic and  solar  data,  auroral  radar  data  from  Churchill,  Canada,  along 
with  instrumentation  and  A-scope  video  tapes  from  PAR.  M&S  Computing 
is  performing  for  BMDSCOM  an  initial  reduction  of  this  data  including  tape 
stripping,  reformatting,  and  data  management  as  well  as  the  graphic  presenta- 
tion of  numerous  auroral  backscatter  maps.  More  advanced  analysis  of  geo- 
magnetic relationships,  morphology,  and  tracking  effects  are  also  being 
performed.  The  volume  and  diversity  of  the  gathered  data  precludes  an 
exhaustive  analysis  of  all  its  aspects.  However,  it  will  be  made  available  to 
the  scientific  community  so  that  additional  analysis  may  be  performed,  if 
desired. 

Section  2 describes  the  auroral  test  plan,  the  1975  auroral  test  and 
the  required  hardware  and  software  changes.  Section  3 describes  the  activ- 
ities and  events  during  the  two  data  gathering  periods.  Section  4 discusses 
the  development  and  features  of  the  aurora  off-line  data  reduction  and 
analysis  software  tools.  This  section  also  explains  the  basic  data  reduction 
process.  Section  5 describes  the  data  stripping  and  processing  which  has 
been  performed.  Included  also  are  examples  of  the  graphic  presentations. 
Section  6 discusses  the  goals  of  analysis  efforts  being  performed  for  BMDSCOM. 
This  section  includes  a description  of  advanced  analysis  software  tools 
being  developed  by  M&S  Computing.  Section  7 outlines  future  auroral  data 
gathering  tasks  to  be  performed  at  PAR  and  explains  new  special  techniques 
being  provided. 

In  this  report,  no  attempt  is  made  to  describe  the  results  of  in-depth 
analysis  of  the  data  presented.  The  analysis  is,  in  fact,  just  beginning. 

The  main  point  is  to  describe  the  study,  the  data  which  has  and  will  be 
collected,  and  the  tools  being  used  to  analyze  that  data. 
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EVALUATION  OF  THE  PAR  AURORAL  EXPERIMENTS 


The  aurora  has  been  of  concern  to  the  PAR  system  designers  and  its 
military  user  since  early  in  its  development.  Several  factors  contribute 
to  this  concern.  Included  are: 

o the  orientation  of  the  phased  array  face  which  is  pointed  almost 
directly  toward  geomagnetic  north, 

o the  UHF  operating  band,  and 

o the  very  high-power  transmitted  and  the  exceptional  sensitivity 
of  the  receiver. 

The  net  effect  of  these  factors  is  that  in  the  presence  of  even  weak  auroral 
activity,  a significant  amount  of  transmitted  energy  will  be  backscattered 
to  the  radar.  For  this  reason,  features  were  designed  into  the  PAR  hard- 
ware and  software  to  minimize  the  effects  of  these  echoes.  It  was  especially 
desired  to  prevent  the  false  reporting  of  offensive  reentry  vehicles.  Later, 
during  the  PAR  installation  and  testing  phase,  experiments  were  performed 
in  the  presence  of  aurora  to  test  and  optimize  these  protective  features. 

The  results  have  proven  satisfactory.  It  was  during  these  first  auroral  tests 
by  Bell  Telephone  Laboratories  in  early  1974,  that  many  of  the  auroral 
blocking  features  were  temporarily  disabled  allowing  a significant  amount  of 
auroral  data  to  be  recorded.  Some  of  this  data  has  been  reduced  and  pre- 
sented in  both  graphic  form  and  as  computer  printout.  However,  it  lacked 
the  necessary  experimental  controls  required  for  detailed  scientific 
analysis.  Moreover,  some  of  the  auroral  blocking  functions  were  still  par- 
tially functional.  The  results  of  these  first  tests  combined  with  the  interest 
of  the  scientific  community  in  the  design  of  new  and  even  more  powerful 
long  range  radars  prompted  further  auroral  studies. 

The  PAR  is  particularly  well-suited  for  auroral  backscatter  studies. 
Among  its  features  are: 

o An  operating  frequency  in  the  450  MHz  band. 

o Phased  array  steering,  providing  rapid  wide  angle  scan 

capabilities. 

o Exceptionally  high-peak  and  average  transmit  powers  allowing 
very  long  range  coverage,  well  over  2000km. 

o Excellent  sensitivity  and  noise  performance. 

o Narrow  receive  beam  <1.  5°. 
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o Excellent  range  resolution < 1.  5 km  coupled  with  very  high 
transmit  pulse  energy,  obtained  by  pulse  compression. 

o Optimum  face  orientation  to  minimize  distortion  and  steering 
loss  in  direction  of  maximum  auroral  backscatter. 

o High  performance  data  processing  and  recording  capabilities. 

o Capability  to  perform  simultaneous  measurements  of  auroral 
backscatter  and  satellite  tracking. 

o Sidelobe  suppression  capabilities  to  minimize  unwanted  side- 
lobe  returns. 

These  features  undoubtedly  made  the  PAR  a potentially  valuable  tool 
for  providing  the  data  to  get  a new  look  at  the  aurora  and  to  bolster  and 
clarify  the  auroral  data  collected  previously  by  the  scientific  community. 

It  offered  the  capability  to  collect  data  in  masses  and  with  a resolution  not 
possible  previously. 

In  the  spring  of  1975,  the  requirement  for  performing  these  studies 
was  provided  in  the  SAFEGUARD  Operational  Experience  Plan.  This  pro- 
gram consisted  of  a wide  variety  of  experiments  which  would  determine  and 
provide  for  the  documentation  of  a number  of  performance,  reliability,  and 
potential  capabilities  of  the  now  fully  operational  SAFEGUARD  System.  In- 
cluded in  this  plan  was  the  general  topic  of  auroral  studies.  M&S  Computing 
was  assigned  the  task  of  determining  the  goals  for  this  study,  developing 
the  test  plan,  specifying  special  hardware  and  software  requirements,  and 
coordinating  the  implementation  of  all  aspects  of  the  plan.  In  addition  to 
those  responsibilities,  M&S  Computing  was  faced  with  the  task  of  coordinating 
the  scheduling  of  aurora  data  collection  with  the  military  user  of  the  SAFE- 
GUARD PAR  Radar.  Difficulty  arose  from  the  low  reliability  associated 
with  the  prediction  of  auroral  events,  a deficiency  compounded  by  the  current 
low  levels  of  sunspot  activity  and  subsequently  infrequent  and  weak  aurora. 

Early  in  1975,  meetings  were  held  at  the  Air  Force  Cambridge  Re- 
search Laboratories  in  an  attempt  to  more  firmly  establish  realizable 
testing  goals.  A further  meeting  was  held  at  the  Air  Force  Global  Weather 
Central  at  Offutt  Air  Force  Base  to  determine  the  auroral  prediction  capa- 
bilities which  could  be  provided.  Also  determined  were  the  types  of  geo- 
magnetic, satellite  photographs,  and  solar  data  which  would  be  available  to 
support  the  radar  data.  As  a result  of  this  meeting,  a channel  of  communica- 
tions was  opened  with  the  USAF  which  would  provide  timely  predictive  data 
to  PAR,  preceding  and  during  the  data  collection  period.  Long  range  pre- 
dictions estimated  a maximum  of  auroral  activity  near  September  21,  1975. 
That  date  was  chosen  as  the  tentative  time  frame  for  the  two  nights  of  data 
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collection  which  would  be  performed.  Further  tie-ins  were  made  with  Finley 
Air  Station,  a UHF  radar  just  south  of  PAR,  and  with  the  Fort  Churchill, 
Manitoba  Research  Range,  operated  by  the  Canadian  National  Research 
Council.  These  facilities  provided  additional  capabilities  for  determining 
the  approach  or  presence  of  auroral  activity;  a critical  capability  since 
PAR  is  so  highly  automated  it  lacks  a plan  position  indicator. 


By  early  April,  1975,  the  basic  goals  of  the  experiment  had  been 
established.  Those  goals  are  stated  briefly  below: 

o Preparation  of  auroral  echo  maps  to  develop  a more  precise 
characterization  of  the  structure  and  time  variation  of  the 
aurora  and  to  provide  a better  understanding  of  the  require- 
ments for  the  auroral  exclusion  region. 


o Analysis  of  the  tracking  errors  which  result  from  signal  propa- 
gation through  an  aurorally  perturbed  ionosphere. 


o Analysis  of  the  sidelobe  blanking  effectiveness  in  an  auroral 
environment. 


o Analysis  of  aurorally  generated  noise  radiation. 

o Correlation  of  radar  auroral  observations  with  forecasting  and 
other  predictive  techniques. 

It  was  determined  that  the  radar  would  scan  as  rapidly  as  possible 
the  region  expected  to  produce  auroral  returns.  This  region  is  determined 
by  the  altitude  of  the  ionospl.  ric  E layer  and  the  location  of  the  radar  and  its 
orientation  with  respect  to  the  earth's  geomagnetic  field.  It  was  decided  that 
the  best  compromise  between  minimizing  the  time  required  to  scan  this  region 
and  still  obtaining  frequent  full  space  observations  would  be  provided  by  a 
sequence  of  nine  consecutive  scans  up  to  20°  elevation  followed  by  one  scan 
of  the  entire  region  observable  by  the  radar.  This  sequence  would  be  initiated 
at  the  onset  of  auroral  activity  and  maintained  until  the  activity  ceased.  It 
was  also  determined  that  a number  of  satellite  calibration  spheres  and  Navy 
Transit  satellites  would  be  tracked  during  the  auroral  event.  Analysis  of 
this  track  data  might  provide  valuable  insight  into  the  effects  of  ionospheric 
disturbances  on  tracking  performance,  especially  in  that  it  would  be  possible 
to  spatially  correlate  these  track  data  with  the  location  of  the  aurora. 


The  development  of  a detailed  test  plan  was  quite  involved  in  that  many 
hardware  and  software  modifications  were  required  to  collect  the  search 
and  track  data.  It  must  be  understood  that  the  PAR  is  a military  radar  con- 
trolled by  a tactical  software  process  designed  to  detect  and  track  only 
satellites  and  offensive  reentry  vehicles  and  to  reject  such  interference 
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as  aurora.  Software  and  hardware  changes,  even  on  a temporary  basis,  are 
thoroughly  screened  and  controlled  and  must  be  exhaustively  tested  to  assure 
that  there  is  no  impairment  to  the  military  capability.  Implementation  of 
these  modifications  is  a lengthly  and  arduous  task.  After  several  meetings 
with  Bell  Telephone  Laboratories,  the  SAFEGUARD  prime  contractor,  a test 
requirements  memorandum  was  prepared  on  May  28,  1975,  stating  the  test 
philosophy  and  basic  capabilities  required  as  well  as  identifying  the  software 
and  hardware  modifications  which  would  be  needed.  The  major  system 
changes  are  given  in  the  following  paragraphs. 

The  tactical  software  generates  a range  blanking  gate  which  suppresses 
video  responses  from  range  intervals  in  the  region  where  auroral  backscatter 
is  expected  to  occur.  This  gate  and  other  constraints  had  to  be  modified 
to  allow  maximum  observation  of  the  aurora.  The  normal  scanning  raster  is 
a complex  irregular  sequence  designed  to  maximize  reentry  target  detection. 
This  pattern  was  changed  to  a regular  non-interlacing  television  type  raster. 
This  facilitates  off-line  data  reduction  and  maximizes  time-space  event 
correlation. 

Target  track  initiation  is  normally  automatic  whenever  a target  meets 
certain  criteria.  The  massive  number  of  auroral  echoes  and  their  diverse 
characteristics  create  a high  probability  of  attempted  automatic  tracking,  an 
undesirable  response  in  this  experiment.  For  this  reason,  the  radar  was 
restricted  from  attempting  any  tracks  except  when  manually  directed.  This 
manual  control  was  reserved  for  the  initiation  of  track  on  the  special  spheres 
and  Transit  satellites.  A special  search  routine,  normally  a part  of  the 
tactical  software,  is  used  to  locate  special  objects.  By  searching  a small 
region,  surrounding  the  desired  target,  it  can  be  detected  and  brought  in 
track.  Software  modifications  were  made  to  extend  the  region  in  which  the 
saellites  could  be  tracked  and  to  keep  the  object  in  a high  data  rate  track 
throughout  the  tracking  sequence. 

The  PAR  data  processor  records  digitized  events  from  the  detected 
return  signal.  During  search,  three  events  are  reported  to  the  processor: 
the  time  at  which  the  video  signal  exceeds  the  noise  threshold,  the  amplitude 
of  each  video  peak  or  zero  slope,  and  the  time  at  which  the  video  falls  below 
the  noise  blanking  threshold.  Each  of  these  events  is  called  a reply.  These 
replies  must  be  sent  to  the  data  processor  where  they  normally  undergo  many 
levels  of  processing  along  with  being  recorded.  In  the  presence  of  aurora, 
there  are  many  video  peaks,  in  fact,  it  has  been  observed  that  sometimes 
replies  are  generated  at  the  maximum  theoretical  radar  response  rate. 

Both  hardware  and  software  techniques  are  built  into  the  tactical  system  to 
prevent  data  overflow.  The  hardware  limits  the  number  of  peaks  which  may 
occur  between  up  and  down  threshold  crossings  while  the  software  limits  the 
total  number  of  replies  which  can  be  generated  during  each  search  cycle. 

The  search  cycle  consists  of  the  transmitting  and  receiving  of  three  pulses  at 


three  separate  angles.  This  software  and  hardware  limiting  was  modified  for 
the  auroral  test  creating  a possibility  of  overloading  the  radar/data  processor 
data  buffer,  the  data  processor,  or  the  recording  equipment.  Overloading 
was  alleviated  by  disabling  one  of  the  three  beams  generated  each  minor  cycle. 
Disabling  one  beam  is  possible  since  there  is  a separate  signal  processor  for 
each  beam  which  permits  the  radar  to  sample  three  directions  simultaneously. 
The  gap  created  by  the  missing  beam  was  closed  by  modifying  the  raster  to 
generate  a different  beamstacking  arrangement  as  shown  in  Figure  2-1. 

Since  the  radar  functions  are  performed  in  timed  sequence  blocks,  the 
number  of  potential  radar  replies  was  cut  by  one -third.  In  addition,  all  PAR 
software  functions  not  needed  for  aurora  data  collection  were  disabled. 

The  PAR  normally  switches  its  operating  frequency  at  regular  periods. 
The  range  of  the  several  frequencies  employed  covers  a broad  band.  A 
technique  was  developed  to  constrain  the  radar  to  operate  at  one  frequency 
at  least  99%  of  the  time. 

The  PAR  receiver  employs  linear  FM  pulse  compression  with  software 
selectable  Doppler  offset  frequency  generators  which,  through  the  mixers, 
compensate  for  target  Doppler.  Doppler  compensation  is  impossible  for  un- 
predictable auroral  measurements  so  modifications  were  made  to  hold  the 
offset  to  zero. 

The  functions  of  the  noise  threshold  control  and  the  IF  attenuator  are 
dynamically  controlled  by  the  computer.  Modifications  were  made  to  allow 
manual  control  of  the  IF  attenuators  and  the  noise  threshold  was  set  to  a fixed 
value.  Similarly,  the  sensitivity  time  control  was  disabled  allowing  full  gain 
from  minimum  to  maximum  range. 

The  last  modification  allowed  the  antenna  pattern  of  the  sidelobe  blank- 
ing system  called  a Q-Channel,  to  be  modified  by  card  reader  input  to  facilitate 
convenient  and  rapid  testing  of  the  sidelobe  blanking  function. 

Guidelines  for  a basic  procedure  for  performing  the  test  were  published 
by  M&S  Computing  on  August  20,  1975.  In  response  to  this  document,  PAR 
Site  personnel  produced  a detailed  test  procedure.  One  of  its  prime  objectives 
was  to  establish  the  procedures  for  initiating  data  collection.  That  is,  terminat- 
ing the  tactical  functions  of  the  PAR,  loading  the  auroral  software,  installing 
the  necessary  hardware  modifications,  and  setting  up  the  special  instrumenta- 
tion and  video  recording  equipment.  One  individual,  a member  of  the  PAR 
Site  Engineering  staff,  was  selected  to  act  as  Test  Director.  The  Test  Director 
was  responsible  for  the  entire  data  collection  process  including  coordination 
with  the  military.  Geomagnetic  indices,  auroral  forecasts,  radar  observations, 
and  visual  observation  data  were  to  be  provided  via  telephone  or  TWX  to  the 
Test  Director.  The  Test  Director  used  those  data  to  decide  whether  to 
interrupt  the  tactical  software  and  implement  auroral  data  collection  and  later. 
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to  decide  when  to  end  the  data  collection.  Once  the  data  collection  process 
was  started,  automatic  recording  was  performed  continuously  until  local  in- 
dicators and  instrumentation  indicated  that  the  auroral  activity  had  subsided 
and  that  the  various  outside  agencies  indicated  a low  probability  of  a re- 
occurrence. 

Special  events  collection  included  the  tracking  of  a number  of  calibra- 
tion spheres  and  Navy  Transit  satellites.  Testing  of  a temporarily  modifed 
sidelobe  blanking  system  or  Q -Channel  would  be  performed  twice  during 
especially  active  periods.  The  auroral  noise  measurement  test  was  planned 
to  consist  of  scanning  at  all  normal  pointing  angles  but  with  the  transmitters 
inhibited.  The  object  was  to  determine  if  there  was  a significant  increase  in 
noise  generated  returns  when  the  radar  was  pointed  toward  the  aurorally 
active  region. 


3. 


GATHERING  THE  AURORAL  DATA 


Site  checkout  of  the  auroral  software  and  hardware  modifications  was 
completed  by  September  13,  1975.  Several  runs  were  made  in  the  presence 
of  moderate  aurora.  The  system  was  functioning  as  planned;  also  the  data 
gathering  procedures  and  alert  scheme  were  being  performed  properly.  The 
tie-ins  to  Global  Weather  Central,  Fort  Churchill,  and  the  Finley  Radar 
Station  were  reliable. 

Global  Weather  Central  had  moved  the  expected  auroral  peak  to  the 
September  15-17,  1975,  timeframe.  The  schedule  change  required  the 
that  the  installation  of  TV  monitors  and  recording  equipment  be  expedited, 
and  the  task  was  completed  by  the  night  of  September  16,  1975.  There  were 
a number  of  video  displays.  One  was  from  a wide  angle  low  light  level  tele- 
vision camera  mounted  on  the  roof  of  the  PAR.  That  camera  was  locked  in 
place  and  oriented  toward  the  geomagnetic  north  pole.  A second  camera 
monitored  an  oscilloscope  which  was  set  up  to  provide  an  A -scope  presenta- 
tion of  the  detected  video.  A third  camera  monitored  a console  Cathode  Ray 
Tube  (CRT)  which  displayed  satellite  ground  tracks  on  a projected  map  of 
North  America.  A fourth  camera  monitored  the  saturation  alarm  counter 
which  indicates  input  signal  levels  above  the  usable  range  of  the  signal 
processor.  That  signal  would  be  used  to  alert  the  Test  Director  that  more 
signal  processor  front-end  attenuation  should  be  added.  The  final  camera 
monitored  the  highly  accurate  system  clock.  Using  the  split  screen  capability, 
this  time  was  patched  onto  a section  of  each  of  the  other  displays,  thus,  pro- 
viding a very  accurate  record  of  time  for  all  recordings. 

It  had  been  observed  at  the  PAR  that  generally  if  auroral  activity  was 
going  to  occur,  it  would  first  appear  between  23:30  and  01:30  UT,  or  18:30 
and  20:30  local  DST.  This  pattern  is  typical  of  the  more  common  diffuse,  quiet 
forms  of  auroral  activity  which  occur  in  periods  of  low  solar  and  geomagnetic 
activity.  Each  evening  the  auroral  software  was  loaded  and  the  A -scope  was 
observed  for  the  unmistakable  signals  telling  of  the  presence  of  aurora. 

3.  1 First  Recorded  Occurrence 


At  midafternoon,  on  September  17,  1975,  Global  Weather  Central  pre- 
dicted a rather  low  probability  of  aurora.  At  22:00  UT,  they  indicated  signifi- 
cant geomagnetic  disturbances  in  Norway  which  increased  the  global  Ap  index. 
At  23:30,  the  auroral  software  was  loaded  into  the  PAR  data  processor,  the 
A-scope  presentation  immediately  showed  the  presence  of  weak  auroral 
echoes.  Unfortunately,  it  was  cloudy  so  no  visual  aurora  was  observed. 

Fort  Churchill,  located  about  1150  km  from  the  PAR  on  a line  almost  directly 
between  the  PAR  and  the  north  geomagnetic  pole,  reported  no  visual  or 
radar  aurora.  The  radar  at  Fort  Churchill,  however,  is  capable  of  only  5kW 
peak  output  power. 


As  observation  continued,  peak  backscattered  signal  strength  as  ob- 
served on  the  A-scope  increased  by  about  11  dB  by  23:54  UT,  then  steadily 
increased  another  5 dB  by  01:00  UT.  The  intensity  of  these  peaks  then  be- 
gan a rapid  increase,  gaining  another  10  dB  by  01:58  UT  and  remained  strong 
until  about  02:30.  When  it  began  to  fall  off,  by  03:15  no  echoes  could  be  ob- 
served. 


At  this  time.  Global  Weather  Central  reported  intense  auroral  activity 
over  Norway.  The  Ap  index  was  up  to  41.  Fort  Churchill  reported  radar 
auroral  activity  and  noted  their  peak  intensity  occurred  about  34  minutes  be- 
fore the  peak  at  the  PAR. 

At  the  PAR,  a few  weak  echoes  were  observed  at  03:50.  These  had 
subsided  by  04:00.  At  04:30,  the  intensity  slowly  increased.  At  05:05,  a 
more  rapid  buildup  was  observed,  but  these  levels  were  still  far  below  the 
earlier  peak.  A steady  buildup  occurred  until  the  second  peak  at  about  07:00 
UT.  The  aurora  then  gradually  decreased  in  intensity  with  frequent  secondary 
peaks  until  10:45  UT  or  5:45  local  when  all  activity  ceased.  Fort  Churchill 
was  at  this  time,  experiencing  a buildup  of  aurora  probably  due  to  the  pole- 
ward  motion  of  the  receding  auroral  oval. 

The  auroral  data  collection  was  terminated  at  10:45  UT.  Two  major 
peaks  had  occurred,  21  different  satellite  tracks  had  been  performed,  and  31 
reels  of  2400  feet  data  tapes  had  been  generated.  The  intensity  of  the  auroral 
echoes  was,  however,  much  lower  than  had  been  hoped  for.  The  auroral 
echoes  failed  to  exceed  the  dynamic  range  of  the  signal  processor.  The 
aurora,  which  was  producing  55-65  dB  s/n  signals  during  the  most  intense 
aurora,  could  not  have  produced  a significant  number  of  detectable  sidelobe 
returns,  so  the  sidelobe  blanking  test  was  not  performed. 

3. 2 Second  Auroral  Occurrence 

The  same  daily  procedure  of  receiving  auroral  predictions  and  early 
evening  radar  observations  was  made  with  no  avail  until  the  night  of  September 
26,  1975.  On  this  evening,  fairly  strong  (+40  dB  s/n)  peaks  were  observed 
as  soon  as  the  auroral  process  was  loaded  at  23:23  UT  (18:23  local),  record- 
ing was  begun  immediately.  A rapid  buildup  in  intensity  was  experienced 
and  there  were  at  least  four  cycles  of  auroral  buildup  and  drop-off  before  the 
major  peak  occurred  near  08:30  UT.  That  peak  yielded  a number  of  echoes 
well  above  75  dB  s/n.  There  was  also  a much  more  pronounced  north/ south 
movement  of  the  southern  edge  of  the  reflecting  region.  At  01:45,  surface 
atmospheric  measurements  were  made  to  be  used  for  lower  troposnheric  re- 
fraction calculations.  At  02:15,  the  PAR  was  experiencing  strong  auroral 
echoes.  Fort  Churchill  reported  only  mild  auroral  activity  and  mild  echoes  with 
a steady  magnetic  field.  Clouds  obscurred  any  visual  aurora  at  the  Canadian 
site.  At  PAR,  the  sky  was  clear  enough  to  see  stars  but  no  visual  aurora  was 
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observed  at  any  time  during  the  night.  At  03:30  UT,  Fort  Churchill  reported 
significant  magnetic  disturbances,  especially  in  the  Z-component.  They  re- 
ported that  these  levels  are  usually  associated  with  a magnetic  storm.  The 
riometer  also  showed  a higher  absorption  level.  At  04:00,  Churchill  re- 
ported similar  conditions  with  moderate  auroral  activity.  PAR  was  ob- 
serving only  moderate  backscatter  levels.  At  04:10,  Global  Weather  Central 
predicted  an  Ap  ^ 30  for  the  next  three  hours.  At  05:00,  they  predicted  an 
Ap  of  approximately  40  for  the  next  six  hours.  The  intensity  of  the  auroral 
returns  was  steadily  increasing,  a broad  fluctuating  plateau  of  strong  returns 
was  experienced  from  06:00  until  07:50  when  the  intensity  began  to  climb 
rapidly.  At  07:15,  the  Finley  Radar  Station  reported  that  it  was  observing 
strong  auroral  returns  from  130  km  to  their  maximum  range  in  a sector  be- 
tween 35°  and  350°  north.  That  meant  that  the  southern  edge  of  the  aurora 
had  moved  behind  the  PAR.  An  intense  peak  was  observed  at  PAR  at  about 
08:30  followed  by  a slight  drop-off  with  another  broad  peak  at  08:49,  followed 
by  a steady  fluctuating  condition  until  10:00,  at  which  time  the  data  collection 
was  terminated.  Thirty-one  data  tapes  were  produced.  Sidelobe  blanking 
tests  were  run  from  07:56:32  until  08:20:05  and  again  from  08:22:10  until 
08:32:10.  The  radar  transmitters  were  inhibited  from  7:48:31  until  7:49:31 
in  an  attempt  to  measure  auroral  noise.  It  should  be  noted,  however,  that 
the  noise  threshold  was  not  lowered  for  this  test,  and  unless  the  auroral 
noise  level  was  very  strong,  probably  very  few  noise  returns  would  be  de- 
tected. A significant  amount  of  satellite  track  data  was  taken  the  night  of 
September  26,  1975. 


4. 


DEVELOPMENT  OF  AURORAL  DATA  REDUCTION  TOOLS 


The  post-processing  of  the  data  from  the  two  data  collection  periods 
is  being  done  at  M&tS  Computing,  Inc.,  in  Huntsville,  Alabama.  Develop- 
ment of  auroral  data  reduction  and  analysis  tools  began  several  weeks  be- 
fore the  data  was  collected.  In  the  planning  stage,  it  was  obvious  that  the 
storage  and  retrieval  of  such  a mass  of  data,  as  produced  by  the  auroral 
test,  would  require  a special  data  management  system.  Tapes  originating 
at  PAR  contain  a large  amount  of  redundant  data  and  the  format  is  not 
suitable  for  production  processing.  Furthermore,  the  data  is  not  organized 
in  sets,  comprising  single  scans  of  the  observed  region  providing  this 
reorganization  was  one  of  the  prime  functions  in  the  initial  data  reduction. 
Another  requirement  was  to  provide  a technique  so  that  the  stripping  and  re- 
reformatting  operation  could  be  made  in  spite  of  occasional  parity  errors 
on  the  original  tapes. 

The  auroral  post -processing  functions  can  be  divided  into  three  major 
blocks:  stripping  of  the  original  data  tapes,  the  mathematical  processing 
functions,  and  the  graphics  display  processor.  A block  diagram  of  these 
functions  is  shown  in  Figure  4-1. 

4.  1 Tape  Stripping  Operation 

The  upper  section  of  Figure  4-1  illustrates  the  overall  functions  of  this 
program.  Auroral  tapes  from  PAR  are  read,  track  data  is  separated  from 
search  data,  redundant  data  is  eliminated,  and  tape  errors  are  corrected. 
Search  data  is  organized  into  scans.  The  stripping  program  numbers  these 
scans,  generates  scan  characterization  data  by  recording  the  start  and  end 
times,  the  number  of  replies  in  the  scan,  and  the  peak  and  the  mean  intensity. 
The  reformatted  search  data  is  then  stored  on  a second  magnetic  tape.  The 
stripping  program  then  stores  the  scan  characterization  data  on  a magnetic 
disk.  Also  recorded  are  reel  identification  number  and  logical  record  number 
of  the  tape  on  which  the  reformatted  search  data  was  stored.  Figure  4-2 
shows  a sample  printout  of  the  data  in  the  search  catalog.  The  catalog  also 
contains  a code  indicating  that  errors  may  have  been  present  on  the  original 
PAR  tapes. 

Satellite  track  stripping  and  cataloging  are  accomplished  in  a similar 
manner.  The  reformatted  track  data  is  stored  on  a separate  set  of  tapes 
than  the  search  data.  An  example  of  the  track  catalog  is  shown  in  Figure  4-3. 

4.  2 Data  Retrieval 

Data  retrieval  is  accomplished  in  the  following  manner.  When  the  user 
has  determined  the  scan(s)  or  track(s)  to  be  processed,  the  scan  (track) 
number  is  input  through  the  desired  mathematical  preprocessor  which  accesses 
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the  search  (track)  catalog.  The  preprocessor  determines  from  the  catalog 
the  basic  information  concerning  the  storage  of  the  pertinent  data.  The 
program  then  requests  the  computer  operator  to  mount  the  proper  tape  reels. 
The  desired  data  is  then  transferred  from  the  tape  to  a disk  drive  for  access 
by  the  computer. 

4.  3 Graphics  Processing 

Graphics  presentation  was  used  whenever  possible  as  a means  to 
output  the  results  of  the  data  reduction  process.  This  technique  was  chosen 
because  of  the  immediate  clarity  and  insight  which  may  be  realized  upon 
examination  of  these  graphical  forms.  Graphical  data  also  is  self -document- 
ing which  is  a real  advantage  when  the  data  is  to  be  distributed. 

One  prime  factor  in  the  choice  of  this  media  was  the  existence  of  a 
sophisticated  Interactive  Graphics  Design  System  (IGDS)  at  M&S  Computing. 
The  extensive  capabilities  of  IGDS  greatly  simplified  the  task  of  providing 
auroral  data  in  a sophisticated  graphical  form. 
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The  IGDS,  which  is  shown  in  Figure  4-4,  allows  the  operator  to 
communicate  with  the  system  by  simple  manual  selections  of  commands  on 
a menu.  Selection  of  these  commands  initiates  the  execution  of  complex 
programs  but  eliminates  the  need  for  time-consuming  card  input.  The 
terminal  keyboard  is  used  when  completely  new  data  is  to  be  processed. 

Key  to  the  operation  of  the  system  are  the  19"  storage  Cathode  Ray 
Tubes  (CRT's).  These  modified  Tektronix  4014  graphics  terminals  have  an 
exceptional  number  of  built-in  capabilities  and  can  be  controlled  to  produce 
complex  plots  at  a very  high  speed.  Dual  screens  are  employed  so  that 
different  plots  can  be  compared  easily.  Data  can  be  quickly  exchanged  be- 
tween the  two  screens.  Permanent  copies  can  be  made  rapidly  using  the  hard 
copy  unit.  If  larger,  multicolored  or  even  higher  quality  copies  are  required, 
output  can  be  directed  to  the  high-speed  CalComp  plotter. 

IGDS  is  an  exceptional  analysis  tool  in  that  it  allows  the  displays  to 
be  instantly  blown  up  or  reduced  to  any  scale.  Another  important  feature  is 
that  the  system  provides  32  overlay  levels  on  each  plot.  Any  combination 
of  these  overlay  levels  can  be  observed  at  one  time.  This  capability  is  used 
to  generate  the  various  altitude  and  azimuth  slices  which  will  be  described 
late  r. 


The  IGDS  Data  Management  System  will  automatically  store  and  catalog 
all  plots  on  mass  storage  devices  so  that  any  item  of  processed  auroral  data 
is  available  for  near  instant  recall. 
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Figure  4-4 


The  auroral  post-processing  effort  has  benefited  greatly  from  the 
availability  of  this  1GDS,  especially  in  the  expense  of  software  development. 
The  Auroral  Graphics  Processor  does  the  basic  data  preparation,  layout, 
and  scaling  as  well  as  supplying  headers  and  scale  information.  From  that 
point,  this  data  is  handled  completely  by  the  IGDS. 

4.  4 Data  Preprocessing 

Preprocessing  involves  the  various,  mathematical,  statistical,  and 
sorting  operations,  which  must  be  accomplished  before  the  data  can  be  pre- 
pared for  graphical  presentation.  There  are  three  types  of  preprocessors, 
as  shown  in  the  middle  section  of  Figure  4-1.  These  are,  from  left  to  right, 
the  Search  Reply  Request  Processor,  the  Activity  and  Event  Profile  (AEP), 
and  the  Track  Analysis  Program. 

4.  5 Activity  and  Event  Profile 

The  AEP  is  a means  for  quickly  determining  the  activity  level  of  the 
aurora  over  the  night  of  the  test.  Three  indicators  of  auroral  activity  are 
used,  number  of  returns  per  scan,  mean  intensity  per  scan,  and  peak  in- 
tensity per  scan.  None  of  these  indicators  completely  describes  the  activity 
level;  so  all  three  are  provided.  Figures  4-5,  4-6,  and  4-7  are  examples  of 
the  three  AEPs  from  September  17  and  18,  1975. 

Each  AEP  plot  contains  one  data  point  for  each  scan.  The  abscissa 
can  be  scaled  to  show  the  entire  night's  activities  or  for  shorter  segments  of 
time,  thereby  allowing  a choice  of  the  time  resolution.  The  AEP  is  gen- 
erated from  the  Search  Catalog  alone  with  no  additional  mathematical  cor- 
rection of  any  kind. 

4. 6 Search  Reply  Processing 

In  planning  for  the  design  of  the  auroral  search  data  processor,  con- 
siderations of  the  mass  of  data  involved  dictated  that  techniques  be  developed 
to  provide  the  researcher  with  means  to  quickly  gain  insight  into  the  overall 
auroral  phenomenon.  The  researcher  could  then  more  readily  exercise  the 
the  judgement  required  to  define  those  periods  in  which  to  concentrate  his 
detailed  analysis. 

One  overview  technique,  AEP,  has  already  been  described.  Another 
means  is  to  generate  auroral  maps  with  the  auroral  echoes  projected  on  the 
earth's  surface  or  some  other  reference  plane.  In  this  way,  the  observer 
can  see  both  the  location  and  relative  intensity  of  the  auroral  backscattering 
from  different  perspectives.  It  was  decided  that  the  metric  system  would  be 
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ACTIVITY  AND  EVENT  PROFILE 


4 TO  1372 


TIME  or  DAY 

ACTIVITY  AND  EVENT  PROFILE 


used  exclusively.  The  different  types  of  auroral  maps  will  be  described  in 
the  next  sections,  but  some  general  explanations  will  be  made  here. 

The  PAR  records  only  certain  portions  of  the  detected  video  signal. 
Recorded  data  are  known  as  replies.  The  three  types  of  events  recorded 
are  noise  blanking  threshold  upcrossings,  video  peaks,  and  noise  threshold 
downcrossings.  The  time  of  each  event  is  also  recorded  as  is  the  log 
amplitude.  Figure  4-8  illustrates  these  events. 

In  all  of  the  auroral  maps  shown  in  this  report,  aurora  is  represented 
by  placing  a character  at  the  projected  location  of  each  video  peak.  For 
each  plot,  the  characters  are  quantized  into  eight  sizes  indicating  eight 
bands  of  log  video  amplitude.  Quantization  intervals  are  determined  using 
the  minimum  and  maximum  reply  on  each  scan.  This  technique  provides 
maximum  resolution  of  echo  strengths  on  each  scan  but  does  not  permit 
direct  comparisons  of  absolute  echo  strengths  of  different  scans. 

Data  for  all  types  of  search  reply  processing  must  undergo  certain 
preliminary  processing  including:  conversion  of  the  radar  processor  bit 
code  into  recognizable  power  level  and  target  range.  The  power  level  term 
is  corrected  for  a number  of  factors  including  transmit  and  receive  antenna 
element  pattern  factor  and  antenna  aperture  area  projected  in  the  direction 
of  the  beam.  A range  to  the  fourth  power  factor  has  been  applied  to  correct 
for  range  attenuation.  This  correction  will  be  changed  to  range/ squared 
when  volume  reflection  properties  are  calculated.  Two  other  factors  are 
also  applied  to  the  signal.  These  are  sensitivity  corrections  determined 
every  12  seconds  during  the  radar  self-calibration  or  Class  A tests.  One 
calibration  factor  corrects  for  changes  in  peak  transmitted  power  level,  the 
other  for  receiver  sensitivity.  Separate  calculations  are  made  for  the  right - 
and  left-hand  beam  in  each  pair  as  separate  signal  processors  are  used  for 
each  one.  Search  processing  may  be  performed  on  either  or  both  beams  so 
that  comparisons  may  be  made. 

Another  function  performed  to  all  search  data  is  coordinate  conversion. 
This  process  converts  the  data  from  the  radar  face  coordinate  systems  to 
N-S,  E-W  projections  or  longitude  and  latitude  with  spherical  altitude. 

One  of  the  capabilities  of  the  basic  search  preprocessor  is  spatial 
data  filtering.  Filtering  may  be  performed  to  reduce  data  point  density  to 
relieve  congestion  on  graphical  auroral  maps,  to  reduce  computer  processing 
time,  and  to  select  data  only  from  special  regions  in  space.  The  ability  to 
select  data  from  special  regions  provides  a powerful  tool  for  scientific 
analysis  of  auroral  data.  The  technique  in  which  this  filtering  approach  was 
implemented  makes  the  overall  processing  of  small  regions  very  efficient  in 
that  it  allows  undesired  data  to  be  bypassed  with  very  little  handling  (a  much 
desired  capability  when  each  scan  contains  as  many  as  27,000  replies). 
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2 - Video  Peak 

3 - Downcrossing 


Figure  4-8 


The  filter  allows  the  user  to  define  a volume  of  space  outside  of  which  no  data 
will  be  processed.  In  addition  to  this  regional  filter,  it  is  possible  to  select 
only  data  from  predetermined  range  intervals  at  any  number  of  selected 
steering  angles.  Efficient  filtering  is  possible  because  all  recorded  data  is 
stored  in  the  same  chronological  order  it  came  from  the  radar.  Radar 
scanning  is  a TV-like  raster  and  replies  are  time  (range)  ordered.  The  spatial 
filter  boundaries  are  converted  to  basic  radar  bit  code  so  that  comparisons 
will  require  no  correction  or  conversion.  Furthermore,  the  sorting  operation 
takes  advantage  of  the  time -space  organizations  of  the  input  data  to  quickly 
eliminate  unwanted  data. 

It  is  often  desirable  to  sort  data  from  different  geomagnetic  off- 
perpendicular  contours.  One  of  the  critical  factors  in  the  auroral  backscatter 
mechanism  is  the  angle  made  between  the  radar  beam  and  the  earth's  magnetic 
field  lines  at  the  reflecting  point  shown  in  Figure  4-9.  Only  when  this  angle 
is  near  0°  (perpendicular)  does  strong  reflection  take  place.  It  is  often  de- 
sirable to  study  the  relationship  of  the  off -perpendicular-to-received  signal 
strength.  The  spatial  filter  can  thus  be  used  to  sort  data  according  to  its 
relationship  with  the  earth's  magnetic  field  if  supplied  with  the  proper  angles 
and  range  intervals.  Information  is  supplied  by  another  program  which  has 
generated  and  stored  this  data  in  tables.  The  program  generates  each  of  the 
discrete  steering  angles  used  while  the  PAR  is  in  the  search  mode.  Then 
using  a 42-term  sperhical  harmonic  solution  of  the  International  Geomagnetic 
Reference  field,  the  program  generates  the  ranges  to  each  point  where  the 
beam  meets  the  desired  aspect  conditions. 

As  a point  of  interest,  the  family  of  all  space  points  where  a radar 
makes  the  same  degree  of  off-perpendicular  angle,  forms  a dome -like 
surface  (open  side  toward  the  earth).  The  set  of  domes  for  all  aspect  angles 
is  stacked  one  inside  the  other  like  mixing  bowls  with  the  0°  off-perpendicular 
surface  on  the  inside.  As  can  be  seen  in  Figure  4-9,  the  aspect  angle 
(90°  -5  ) can  be  both  positive  and  negative.  Because  of  their  location,  most 
radars  cannot  make  negative  aspect  angles  because  these  contours  are  stacked 
inside  the  0°  surface  and  are  observable  only  at  very  low  elevations.  The 
earth's  curvature,  then,  usually  blocks  observation  of  these  regions.  In 
fact,  the  PAR  can  never  make  less  than  +1.0°  off-perpendicular  in  the 
altitude  region  where  aurora  occurs. 

The  tables  generated  by  the  geomagnetic  model  can  be  directly 
accessed  by  the  spatial  filter  to  select  data  which  corresponds  to  each  off- 
perpendicular  contour.  The  data  can  then  be  used  for  various  analytical 
processing. 
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Figure  4-9 


Several  specific  search  data  processing  capabilities  have  been  planned, 
developed,  and  are  being  used.  Some  capabilities  are  still  in  development  and 
a few  are  still  in  the  preliminary  design  stage.  The  next  sections  describe 
those  programs  which  are  either  fully  operational  or  are  near  completion. 

4.  7 Auroral  Maps 

Three  types  of  auroral  maps  are  provided:  top-down  views,  profile, 
and  sine  space  views.  These  maps  are  presented  on  M&S  Computing's  IGDS. 

All  auroral  plot  data  is  stored  on  magnetic  tape  for  rapid  recall.  Selected 
plots  are  maintained  on  disk  files  for  immediate  access. 

4.7.1  Top-Down  Views 

Top-down  views  provide  a picture  of  auroral  reflections  as  seen  look- 
ing down  from  above  the  aurora  toward  the  center  of  the  earth  with  the  reflec- 
tion points  projected  onto  the  earth's  surface.  The  surface  is  represented  as 
either  a map  of  North  America  or  a rectangular  projection  in  N-S,  E-W 
coordinates.  The  auroral  reflections  are  represented  by  a character  of  size 
proportioned  to  the  logarithm  of  the  corrected  echo  intensity.  Figure  4-10  is 
an  example  of  a top-down  view  with  the  N-S,  E-W  projection.  Lines  extend 
outward  from  the  point  where  the  PAR  is  located,  the  dashed-line  points 
toward  the  geomagnetic  north  pole,  the  vertical  solid  line  is  true -north,  and 
the  solid  lines  are  separated  by  15°.  Top-downs  are  generated  in  sets  of  10 
on  each  plot  corresponding  to  different  slices  in  altitude.  Once  generated, 
these  slices  are  on  separate  levels  and  can  be  observed  singularly  or  in  any 
combination.  Figure  4-10  is  a composite  of  all  altitudes  from  70  to  170  km. 

Figure  4-11  shows  an  auroral  top-down  projected  on  a map  of  North 
America.  The  map,  data,  scales,  coordinates,  geographic  locations,  etc.  , 
are  each  on  separate  levels  so  they  can  be  put  together  in  any  combination. 

4.7.2  Profile  Views 

The  profile  views  involve  a projection  where  the  aurora  is  observed 
form  an  angle  perpendicular  to  the  bisector  of  the  azimuth  slice(s).  The  data 
points  on  the  projection  are  shown  with  their  true  range  because  the  data 
points  in  the  slice  are  folded  into  the  projection  plane  at  the  bisector. 

Figure  4-12  is  a composite  of  all  azimuths  from  -30°  to  +60°. 

4.  7.  3 Sine  Space 

Sine  space  is  a coordinate  system  used  to  represent  radar  pointing  angles 
using  direction  cosines.  Sine  a is  an  azimuth-like  coordinate  and  sine0  is  an 
elevation-like  coordinate.  Those  familiar  with  the  technique  will  find  the  plots 
Informative.  Three  types  of  sine  space  plots  are  provided  with  each  level 
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TOP-DOWN  VIEW  OF  AURORA  PROJECTED  ON  MAP  OF  NORTH  AMERICA 
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displaying  a different  indicator  of  auroral  intensity.  A character  is  located 
at  each  beam  position  where  at  least  one  auroral  echo  was  detected.  The 
size  of  the  character  is  proportional  to  either  the  number  of  peaks  measured, 
the  average  of  all  peaks  from  that  angle,  or  the  maximum  peak.  Figure  4-13 
shows  a typical  sine  space  plot. 

In  the  upper  right-hand  corner  of  all  maps  is  a legend  telling  whether 
the  right-  or  left-hand  beam  is  included,  the  scan  number,  and  the  Julian 
data  (to  one  second  resolution)  of  the  start  and  end  time  of  this  scan. 

4.7.4  Histograms 

The  histogram  preprocessor  does  not  generate  maps,  but  will  be  dis- 
cussed here.  The  histogram  program  generates  frequency  of  occurrence 
statistics  for  various  data.  Almost  any  set  of  parameters  could  be  processed 
by  this  program.  Figure  4-14  shows  a histogram  of  the  corrected  signal-to- 
noise  ratio  of  all  replies  in  scan  190. 

4.  8 Track  Data  Processing 

There  are  a number  of  programs  included  in  the  track  data  processing 
package.  Some  require  only  the  raw  unsmoothed  radar  track  data,  others 
require  the  input  of  precise  satellite  trajectory  information,  known  as 
ephemeris  data.  Two  primary  factors  have  been  chosen  in  the  study  of  aurora 
and  tracking  interactions.  These  are  scintillation  and  refraction.  Track 
processing  software  which  is  fully  developed,  is  designed  for  the  analysis  of 
scintillation  effects.  There  are  plans,  however,  to  extend  the  analysis  to 
include  refraction  studies. 

Scintillation  appears  as  rapid  fluctuation  in  signal  power  and  apparent 
target  position  which  can  result  from  perturbations  of  the  radar-to-satellite 
propagation  path.  It  is  generally  accepted  that  these  perturbations  are  caused 
by  irregularities  in  the  ionospheric  F region.  It  is  further  apparent  that  the 
characteristics  of  these  irregularities  are  related  to  the  magnetic  disturbances 
associated  with  auroral  activity. 

Scintillation  may  be  studied  in  one  way  by  analyzing  the  point-by-point 
fluctuations  in  signal  strength  and  position.  For  this  study,  absolute  target 
position  is  not  important  and  radar  biases  and  refraction  errors  are  negligible 
because  these  factors  change  slowly  and  smoothly. 

Two  programs  are  presently  available  for  the  analysis  of  ionospheric 
scintillation.  For  these  programs,  equations  have  been  defined  which  process 
the  point -by-point  raw  unfiltered  radar  track  data  and  produce  index  terms 
which  are  related  to  the  fluctuations  in  the  received  signals. 
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Fluctuations  in  sum  beam  amplitude  are  characterized  by  the  dimen- 
sionless Scintillation  Index  (SI),  given  by: 

Sip  = Pmax  ~ Pm  in  4 ^ 

Pmax  + pmin 

Pmax  an<I  I^min  are  t^ie  maximum  and  minimum  power  in  watts  measured 
during  given  time  intervals.  The  radar  tracks  at  a constant  pulse  rate  so 
this  interval  may  be  set  to  include,  for  example,  five  points.  A sliding 
window  is  used  so  that  the  SI  can  be  calculated  with  each  new  track  data  point 
by  including  the  last  four  points  which  were  used  to  calculate  the  previous 
SI. 


The  effect  of  the  Sip  is  to  express  the  signal  fluctuation  relative  to  the 
average  signal  strength.  An  example  is  shown  in  Figure  4-15.  Here,  the 
ordinate  is  the  Sip  and  the  abscissa  is  the  time  into  track  referenced  to  the 
satellite  acquisition  time.  The  amplitude  SI  is  most  meaningful  for  satellites 
with  a uniform  cross  section  such  as  spheres.  The  example  is  for  satellite 
number  4963,  a sphere  with  a .26m  diameter  (note  the  sudden  bursts  of 
scintillation).  These  bursts  have  been  observed  with  many  satellites  where 
the  radar  propagation  path  passed  above  the  auroral  region.  Another  program, 
described  later,  can  be  used  to  correlate  this  scintillation  with  the  propagation 
path. 


Fluctuations  in  the  radar  angle  measurements  cannot  be  characterized 
in  the  same  way  as  power  and  range  because  angle  errors  are  absolute  and 
their  degree  of  degradation  is  generally  not  related  to  the  magnitude  of  the 
steering  angle.  Angular  fluctuations  are,  thus,  expressed  in  terms  of  their 
Perturbation  Index  (PI),  given  by: 


PI 


'max 


®min 


4.2 


The  PI  can  be  calculated  for  different  types  of  angles,  such  as  the  off-bore- 
sight  angular  displacements  a sin  a , A sin  0 , or  for  the  actual  azimuth  and 
elevation  measurement.  The  PI  is  calculated  using  a sliding  window  just  as 
for  the  SI.  Figures  4-16  and  4-17  are  a plot  of  PI  for  satellite  2759.  The 
ordinate  of  these  plots  shows  the  PI  for  Asincr  and  Asin/3  , respectively. 
Figures  4-18  and  4-19  give  the  PI  for  azimuth  and  elevation.  The  PRC  = 
SIPT05  indicates  that  a five  data  point  or  one  second  window  was  used. 


A perturbation  index  PIR  is  also  calculated  for  fluctuations  in  range. 
The  PIR  is  given  by: 
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PERTURBATION  INDEX  FOR  SINcr  ANGULAR  ERROR 
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PERTURBATION  INDEX  FOR  SIN  P ANGULAR  ERROR 
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PERTURBATION  INDEX  FOR  THE  TARGET  AZIMUTH 
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PERTURBATION  INDEX  FOR  THE  TARGET  ELEVATION 
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Another  set  of  programs  is  available  which  compares  raw  radar  data 
to  either  Kalman  smoothed  or  to  precise  ephemeris  data.  The  Track  Analysis 
Program  calculates  residuals  at  each  track  data  point.  Residuals  can  be 
determined  using  either  the  Kalman  smoothed  radar  data  or  the  ephemeris 
data.  Residuals  not  only  show  scintillation,  but  the  effects  of  tropospheric  and 
ionospheric  refraction,  biases  in  the  radar  equipment,  and  bias  errors  in- 
troduced by  coordinate  conversions  due  to  inaccuracies  in  modeling  the  location 
and  orientation  of  the  radar  face. 

In  this  initial  track  data  analysis,  no  attempt  was  made  to  correct  for 
refraction  and  biases.  In  the  best  of  circumstances,  these  effects  can  never  be 
completely  removed  from  the  data.  There  are  a number  of  conflicting  opinions 
and  just  as  many  sophisticated  schemes  for  obtaining  the  best  fit  to  the 
ephemeris  data.  Most  of  these,  however,  require  that  assumptions  be  made 
about  environmental  conditions  some  of  which  are  invalid  during  auroral  activity. 
This  study,  however,  avoids  these  problems  because  these  biases  and  refrac- 
tion errors  can  be  easily  separated  from  the  type  of  scintillation  being  studied. 
This  type  of  scintillation  is  characterized  by  a sudden  increase  in  the  level  of 
signal  fluctuations,  whereas  biases  and  refraction  errors  add  either  a constant  or 
very  slowly  changing  shift  to  the  measured  signal.  The  time  scale  involved  with 
perturbations  due  to  scintillation  and  those  for  refraction  and  bias  errors  are 
separated  by  several  orders  of  magnitude.  If  these  expectations  should  prove 
erroneous,  the  desired  corrections  can  be  readily  applied. 

Two  types  of  programs  are  available  for  performing  scintillation 
analysis  using  residuals.  One  generates  a plot  of  residuals  versus  time  in 
track.  The  residual  can  be  in  sine  space  and  range  or  in  the  radar  face  X,  Y, 

Z coordinate  system.  Another  capability  exists  to  calculate  the  ratio  of  each 
residual  to  the  vector  magnitude  of  the  corresponding  X,  Y,  or  Z component  of 
its  location  at  each  point  in  time. 

Figure  4-20  shows  the  ephemeris  trajectory  of  satellite  3133  in  the  radar 
face  X,  Y,  Z coordinate  system.  Figure  4-21  shows  the  Xpp  track  residuals  for 
this  satellite.  No  attempt  has  been  made  to  correct  for  refraction.  Figure  4-22 
shows  the  same  residuals  with  the  Sip  overlaid. 

The  Track  Analysis  Program  also  generates  geodetic  ground  track  data 
for  each  satellite.  This  information,  which  is  given  in  longitude  and  latitude, 
can  be  used  to  plot  the  path  of  the  satellite  projected  onto  a map  of  North 
America.  In  addition  to  the  satellite  ground  track,  this  program  also  gener- 
ates a ground  track  of  the  entry  and  exit  points  of  the  radar  beam  for  both  the 
ionospheric  E and  F layers.  The  bands  formed  when  these  points  are  plotted 
will  be  used  to  correlate  the  regions  where  scintillation  occurred  with  the  loca- 
tion of  the  0"rora  and  the  auroral  oval  developed  by  Feldstein.  An  example  is 
shown  in  Figure  4-23. 
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The  original  PAR  auroral  tapes  were  copied  on-site  and  mailed  to 
Huntsville.  Data  reduction  was  started  immediately  upon  their  arrival. 

5.  1 Examples  of  Processing  Data  from  September  17  and  18,  1975 

This  section  includes  a large  set  of  self-explanatory  plots.  When 
examining  them,  it  may  be  interesting  to  refer  again  to  the  AEP  in  Section  4.  4. 

A considerable  amount  of  additional  data  of  the  type  shown  has  been 
processed  and  is  available  upon  request.  All  altitude  and  azimuth  slices  indi- 
cated on  the  plots  are  also  available. 

Some  additional  information  about  the  operation  of  the  radar  and  the 
nature  of  the  auroral  backscatter  mechanism  will  be  helpful  in  interpreting 
these  plots.  The  auroral  plots  exhibit  certain  characteristic  shapes  especially 
noticeable  in  the  profile  views.  Some  of  the  outer  limits  of  the  auroral  re- 
flecting mass  are  truly  defined  by  the  structure  of  the  aurora  itself,  but  some 
of  these  boundaries  result  from  radar  scan  limits  or  magnetic  aspect. 

Figure  5-1  will  be  used  to  illustrate  the  mechanism  responsible  for  these 
boundaries. 

One  boundary  is  established  by  the  lowest  elevation  scanned  by  the 
radar.  This  forms  the  maximum  observable  range  boundary.  This  lower 
scan  limit  also  cuts  off  the  outermost  range  extent  in  the  top-down  views. 

The  upper  altitude  bound  of  auroral  backscatter  is  limited  by  magnetic  aspect 
which,  in  the  case  of  the  PAR,  is  no  better  than  5.  5°  at  150  km  altitude  and 
increases  by  roughly  . 5°  per  10  km  altitude  in  the  upper  region  of  observed 
radar  aurora.  The  edge  of  the  reflecting  aurora  facing  the  radar  is  bounded 
by  both  aspect  and  the  lower  edge  of  the  aurora  itself.  The  minimum  radar 
range  used  in  auroral  mapping  very  seldom  affects  the  data  in  any  way. 

In  the  top-down  plots,  the  characteristic  shape  of  the  auroral  reflecting 
region  is  determined  primarily  by  the  lower  elevation  limit  and  aspect.  The 
far  side  of  the  aurora  is  limited  by  the  elevation  limits  as  is  evident  from 
Figure  5-2,  the  near  side  is  limited  by  magnetic  aspect. 

Establishing  these  bounds  is  only  a means  for  interpreting  certain  limit- 
ing criteria  on  the  plots.  In  many  cases,  the  edges  of  the  aurora  occur  in  re- 
gions where  no  limiting  criteria  should  be  of  any  significant  effect.  The  plot 
in  Figure  5-2  is  an  example.  In  this  case,  the  backside  of  the  aurora  is  scan- 
limited  but  at  near  +8°  azimuth  the  near  side  of  the  aurora  is  only  a few 
degrees  off-  X,  At  near  5°  off -perpendicular  the  upper  altitude  region  is  still 
not  severely  affected  by  aspect.  This  example  suggests  that  the  near  and  topside 
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actually  represent  the  physical  limits  of  the  ionospheric  disturbance  that 
causes  the  reflections. 

Figure  5-3  illustrates  the  general  affect  of  aspect.  This  set  of  off- 
perpendicular  contours  approximates  those  that  would  be  expected  for  PAR 
at  a constant  100  km  altitude.  In  the  top-down  views,  the  forward  swept  near- 
side of  the  auroral  reflections  is  a result  of  changes  in  aspect  shown  in  the 
figure.  Any  banding  effects  or  voids  in  the  aurora  normally  truly  represent 
the  true  physical  phenomenon  in  that  they  appear  to  be  unperturbed  regions. 

Figures  5-4  through  5-85  were  generated  from  data  taken  September  17 
and  18,  1975. 

5.  2 Examples  of  Processed  Data  from  September  26  and  27,  1975 

Figures  5-86  through  5-180  are  plots  generated  from  data  taken  during 
the  experiments  on  September  26  and  27,  1975.  The  data  represent  samples 
taken,  in  most  cases,  at  approximately  15-minute  intervals  throughout  the 
data  collection  period.  ,■ 
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6. 


GOALS  OF  THE  ANALYSIS 


The  initial  goals  of  the  study  included: 

o Accurately  defining  the  region  in  space  where  auroral  reflections 
could  affect  the  operation  of  the  PAR  radar. 

o Determining  the  signal  strength  distribution  and  number  of  re- 
plies which  could  be  expected  from  various  regions. 

o Evaluating  the  reliability  of  auroral  predictions  from  various 
outside  sources. 

o Analyzing  the  effects  of  auroral  activity  on  tracking  performance. 

o Characterizing  certain  features  of  the  morphology  of  entire 
auroral  substorms. 

o Determining  if  the  aurora  generated  significant  noise  in  the  PAR 
operating  band. 

o Evaluating  sidelobe  blanking  performance  in  an  auroral  environ- 
ment. 

As  stated  in  the  introduction,  these  analyses  are  incomplete  at  this 
time,  however,  some  very  interesting  data  has  been  presented  here  and  a 
number  of  very  useful  analysis  tools  have  been  described.  Within  a short 
time,  the  results  of  some  of  these  study  areas  should  be  available  and  will  be 
published  at  that  time. 

The  goals  for  analysis  have  expanded  somewhat.  One  of  the  new 
efforts  involves  determining  the  geomagnetic  aspect  dependence  of  the  auroral 
reflections.  This  task  will  be  accomplished  using  the  Geomagnetic  Model  and 
Spatial  Filter  described  in  Section  4.  The  basic  analytical  process  will  be 
similar  to  that  used  by  previous  researchers  except  that  there  will  be  a great 
deal  more  normalization  of  radar  parameters  to  eliminate  the  effects  of  range, 
steering,  and  beam  divergence.  This  normalization  will  be  accomplished  by 
representing  the  auroral  intensity  by  its  volume  reflecting  properties  using  a 
parameter  which  will  be  called  Reflectivity.  This  technique  is  presently  in 
the  final  stages  of  development  at  M&S  Computing. 

Another  new  study  effort  is  directed  toward  determining  the  effect  of 
altitude  on  the  intensity  of  auroral  echoes.  This  will  be  done  by  determining 
the  distribution  of  auroral  reflectivities  from  data  points  at  various  altitudes 
but  on  the  same  geomagnetic  aspect  contour. 
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